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Abstract. We consider laboratory experiments that can detect stable, neutral strongly 
interacting massive particles (SIMPs). We explore the SIMP annihilation cross section 
from its minimum value (restricted by cosmological bounds) to the barn range, and 
vary the mass values from a GeV to a TeV. We also consider the prospects and problems 
of detecting such particles at the Tevatron. 

1 Introduction 

Strongly Interacting Massive Particle (SIMPS), by which we will always mean 
. neutral, stable SIMPs, are of current interest for at least three reasons: 

o 

. • They could be a dark matter constituent as suggested some time ago by 

| Dover, Gaisser and Steigman and by Wolfram Starkman et al.,|| 

r-j ■ show SIMPs would be restricted to rather narrow mass ranges if they were 

P> { to exhaust Q = 1. We will not make this assumption and will consider SIMPs 

outside the regions allowed by the analysis of ref. ||. 
| • It is possible that the lightest SUSY particle (LSP) is strongly interacting and 

. hence, if R-parity is conserved, would form a colorless SIMP. Possibilities, 

such as a gg bound state are discussed in ref. 
• An explanation of the ultra high energy cosmic ray events (UHECRs) pro- 
rN \ posed by Farrar, Kolb and co-workers § is that they are due to interactions 

of SIMPs with a mass below 50 GeV and a cross section for interactions with 
nucleons on the order of a few millibarns or more. 



This summary will review two laboratory experiments that might detect SIMPs. 
More detail can be found in ref. || and the paper on which it is based, ref 0. 

In Section 2 we consider the possibility of finding SIMPs bound in ordinary 
nuclei by searching for anomalously heavy isotopes of higlr-Z nuclei. It is a plea- 
sure to note that the accelerator mass spectrometry (AMS) group at Purdue is 
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in the process of performing the experiment^ suggested in ref. |6j. In Section 3, 
we address the extent to which production and detection of SIMP-anti-SIMP 
(SS) pairs might be performed at the Tevatron. 

Our results, in brief, are that the AMS experiment should be sensitive to 
SIMPs over a wide range of parameter space: (<Jsn, Ms), where Ms is the SIMP 
mass and <jsn is its cross section for scattering off nucleons. The Tevatron, on 
the other hand is likely only to produce and to detect SIMPs in a much more 
restricted range, but one that includes much of the mass range for which the 
SIMP could be the UHECR explanation. It would be only fitting, since much of 
the work on that possibility || was done at Fermilab, if SIMPs were to be de- 
tected at Fermilab and we encourage those with influence in the collaborations 
to explore vigorously that possibility. Finally, we note that we proceed with- 
out committing to a specific SIMP model. We parameterize the experimental 
predictions in terms of the two parameters <Jsn and Ms- 



2 SIMPs in Nuclei 

We know a fair amount about SIMP binding in nuclei from the phenomenology 
of hyper- fragments. See, for example, Povh || for a readable review. Based on 
that experience, we can write for the binding B of the SIMP in a nucleus A the 
relation: 

B= \V SN \-Tr 2 /(2^R 2 ) , (1) 

where fi is the reduced mass of the S-A system, R is the radius of the nucleus A, 
and V is the S-N potential averaged over the volume of the nucleus X. We expect 
the low energy potential, Vsn, to be always attractive. This is true if exchange 
of vacuum quantum numbers dominates. We assume this to be the case, and 
have not found a model to the contrary. Under this assumption, the SIMP can 
be bound in a nucleus for which /i and R 2 ~ A 2 / 3 are large enough to make the 
kinetic energy less than the (average) magnitude of the attractive potential. 

From equation (1) we see that the best chance of finding SIMPs is to search 
in high Z (large) nuclei which minimize the kinetic energy term. Capture by light 
elements at the time of cosmic nucleosynthesis has been studied in ref. || . Atomic 
Mass Spectrometer (AMS) searches to date are reviewed in the careful study of 
Hemmick et al., [ flO| where one learns the somewhat surprising fact that previous 
searches have only been conducted up to sodium (Z = 11, A = 23). This makes 
the current Purdue AMS experiment particularly exciting. They are looking in 
gold (Z = 79, A = 197). 

How big is the potential Vxjv? We take this as a parameter, but we can put 
an approximate LOWER bound on it from the requirement that primordial S 
and S, left over from the early universe, not overdose the universe so that it 
couldn't have continued expanding until today (early 2000). The classic book 

1 We are grateful to Professor Ephraim Fischbach for keeping us informed as to 
progress on this experiment. 
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of Kolb and Turner [fTl]] tells us that the number density of primordial SIMPs 
behaves as 



n s ~ (M s a S g) 



(2) 



Equation (2) says that too small an annihilation cross section means too 
many SIMPs will be left over from the early universe, and Kolb and Turner 
collect together the numerical recipes for computing how small is too small. We 
still need, however, to relate the annihilation cross section, a S g to the SIMP- 
nucleon cross section, usn and to the S — N potential in Equation (1). We make 
the simple ansatz 



where (3 should be on the order of one. Note that Vsn goes as /3 1 / 4 so that our 
results for binding will not be highly dependent on the precision of Equation (3) . 

Now that we know, for each point in the Ms, usn parameter space, the 
primordial S abundance and the binding energy in nuclei, we are almost ready 
to compute for our friends at Purdue, the abundance of anomalous gold-gold 
with a SIMP bound in it. First, however, we need a scenario for how the SIMPs 
get bound into the gold. Our picture is as follows: 

• We assume that the ratio of SIMPs to protons in the galaxy is the same as 
the cosmic ratio, but that most of the SIMPs are in the galactic halo (i.e., 
that their density distribution is p ~ R~ 2 , where R is the distance from the 
galactic center), not in stars. We can then calculate the SIMP flux on the 
Earth, since we know that the Earth is traveling through the galaxy with a 
velocity of about 200km /s which not too different from the galactic virial 
velocity. 

• We assume that when the SIMP hits the Earth, it is slowed by scattering 
with all nucleons and nuclei at a rate determined by cjsn, but can only be 
captured by a nucleus that is large enough. 

• Gold must compete, for SIMP capture, with the most abundant nuclei large 
enough to bind the SIMP. Our comparative estimates use, as the most abun- 
dant elements:, aluminum (A = 27). barium (.4 = 137), and lead (A = 206). 

Our procedure is then as follows: 

• We chose values for Ms and ctsn and then determine whether, for that point 
in parameter space, there is binding in gold. 

• Assuming that there is binding, we then determine (a) the mean free path 
in Earth from the galactic virial velocity and ctsn, and (b) which of the 3 
elements above is gold's chief competitor for SIMP capture. 



Vsn — Vnn (vsn/vnn) 1 



(3) 



a SN = P&NN&ss 



(4) 
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• From the ratio of the abundance of gold to its chief competitor, the mean 
free path, and the average density of Earth, we then compute the chance 
of a particular gold nucleus within a mean free path to capture an incident 
SIMP. Multiplying by the flux (see above) of SIMPs and the time for which 
the sample being put in the AMS target has been exposed gives us the 
fraction of gold nuclei in the sample that should have a SIMP if they exist 
at that point in parameter space. 

Finally, we assume^ that the exposure time is 10 million years because there 
are regions that are geologically inactive over such periods and have had for 
example "placer" gold in the beds of streams for a longer period than that. 

The results are shown in the table. It gives log 10 of the ratio of normal to 
anomalous gold nuclei. The dashes indicate parameter values for which there is 
either no binding in gold or overclosure of the universe. One sees that smaller 
values of asN give larger ratios of anomalous to normal gold. This is because 
smaller values imply that only lead has a nucleus large enough to compete with 
gold for SIMP capture and because the smaller cross section means more pri- 
mordial abundance. The important thing to take away from hours of table study 
is the fact that the relative abundance entries are all considerably higher (for 
anomalous to normal) than the limits of 10 -20 that have been set in AMS work 
on some of the light elements. This provides reason to expect that, if the SIMPs 
are there, the Purdue AMS people will find them. 

3 SIMPs at Fermilab 

Next we consider SS production at the Tevatron. Since we are talking neutral 
SIMPs, we expect little or no signal in the central tracker and in the electro- 
magnetic calorimeter. However, in the hadron calorimeter, we expect to detect 
SIMP signals if <tsn is large enough. The detection of SIMPs is possible if one 
triggers on two relatively back-to-back hadron calorimeter showers, accompanied 
by little else. We will use 10 GeV for the minimum size showers for which such 
triggering might be done. Our task now is to determine: 

• For what values of {Ms, erg} will the SIMP interact in the steel plates of the 
hadron calorimeter? 

• For what values of these parameters will we get calorimeter showers greater 
than 10 GeV or more? 

• Can one recognize a SIMP shower if one sees one? 

• How many such events should we expect? 

First we look at the region of parameter space for which there will be interaction. 
The minimum annihilation cross section permitted by the cosmology argument 

2 We appreciate conversations with Professor E. T. Herrin on searching for old exposed 
gold, and we note that SMU geologist, Dr. Douglas Oliver, has secured such samples 
for the Purdue experiment. 
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Table 1. Mx (vertical) is in units of GeV, and axN (horizontal) is in units of mb. 
Table entries are log 10 (l//), and the — indicates those cases for which X does not bind 
at all. 
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Fig. 1. Contours for energy loss as a function of {Mx , Ex}- The contours displayed 
are in steps of 10 GeV. 

is ~ 3 x 10~ 13 barns, which corresponds through Equation (3) to about a mi- 
crobarn for the S-N cross section. SIMPs with such small cross sections won't 
shower in 1 meter of steel, but for a higher cross section of a few millibarns, we 
would expect 10 or more interactions with the 10 27 nucl eons / cm 2 in the 1 meter. 

To estimate the energy we expect in a shower resulting from a SIMP interac- 
tion in the steel plates of a hadron calorimeter we use a cosmic ray rule of thumb 
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kindly provided by G. Yodr|^ who says that, in a high energy strong interaction, 
about half the center of mass energy goes into inelasticity. In the figure, we give 
the (laboratory) energy released in the calorimeter as we vary the mass and en- 
ergy of the SIMP; the straight lines are constant shower energies. One sees that 
the bigger the SIMP lab energy, the greater a SIMP mass will result in a given 
shower energy. 

Consider now the question of whether we would recognize a SIMP shower if 
we saw one. The background for SIMP showers would likely be neutron showers 
and K decays. The distinguishing feature would be shower opening angle. A 
pion moving transverse in the cm. system would have a lab angle given by 
tan# = I/7. Comparing the angle for a SIMP with that from a neutron of 
the same energy, the SIMP shower should be wider by roughly the ratio of the 
masses. 

Finally, we turn to the number of SIMP pairs the Tevatron might produce. 
We scale the (known) production rate of jets by the ratio of the S-N cross 
section to that of Meson-N, which we take to be on the order of 30 millibarns. 
So long as the SIMP energy is a few times its mass, we don't worry about phase 
space suppression. We assume conservatively a cross section of about 3pb for any 
one parton in the region E > 200GeV. This implies about 6000 events at the 
Tevatron Run II. The estimate of §] is that the Nucleon-UHECR cross section 
needs to be over a tenth the Meson-Nucleon cross section, so we estimate over 
600 events in the Tevatron Run II if SIMPs are the explanation for the UHECR 
events. 

4 Summary 

From the Table we see that there is SIMP binding in gold for MgUsN > 
5mbGeV 2 , and that AMS experiments sensitive to one part in 10 20 can de- 
tect the existence of SIMPs of mass less than a TeV, while the region of interest 
for explaining UHECRs can be explored with a sensitivity of one part in 10 16 or 
less. Looking for SIMPs at the Tevatron is more difficult, but over half the region 
of interest for explaining UHECRs could be searched in the upcoming Run II by 
looking for (wide) back to back jets with no signal in the central tracker or EM 
calorimeter. 

We thank D. Berley, K. Dc, D. Dicus, M.A. Doncheski, R. Ellsworth, G. Far- 
rar, T. Ferbel, E. Fischbach, E.T. Herrin, G. Landsberg, D. Oliver, D. Rosen- 
baum, R. Scalise, and G. Yodh. The work of RNM has been supported by the 
National Science Foundation grant under no. PHY-9802551 . The work of Olness, 
Stroynowski, and Teplitz is supported by DOE. 
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